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We study many interacting Brownian particles under a tilted periodic potential. We numerically 
measure the linear response coefficient of the density field by applying a slowly varying potential 
transversal to the tilted direction. In equilibrium cases, the linear response coefficient is related to the 
intensity of density fluctuations in a universal manner, which is called a fluctuation-response relation. 
We then report numerical evidence that this relation holds even in non-equilibrium cases. This result 
suggests that Einstein's formula on density fluctuations can be extended to driven diffusive systems 
when the slowly varying potential is applied in a direction transversal to the driving force. 
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I. INTRODUCTION 

The construction of statistical mechanics applicable to 
non-equilibrium systems has been attempted for a long 
time. However, even in non-equilibrium steady states, 
no general framework is known, except for the linear 
response theory For example, a useful form of the 
steady state distribution function has not been deter- 
mined in systems far from equilibrium. 

In order to consider a strategy for investigating the sta- 
tistical mechanics of non-equilibrium steady states, let us 
recall the justification of equilibrium statistical mechan- 
ics. The equilibrium statistical mechanics leads to Ein- 
stein's formula on macroscopic fluctuations, which relates 
a large deviation functional of macroscopic fluctuations 
with a thermodynamic function. Since the validity of this 
formula can be checked by measuring the fluctuations of 
macroscopic quantities, we expect that its examination 
in non-equilibrium steady states might provide a hint to 
develop non-equilibrium statistical mechanics. 

With regard to the large deviation functionals in non- 
equilibrium steady states, we review a few works. First, 
a large deviation functional of density fluctuations is ex- 
actly derived for non-equilibrium lattice gases [4, 5]. This 
functional is nonlocal in contrast to equilibrium cases, 
and the local part takes the same form as the equilib- 
rium form. The latter property seems to be specific to 
the model, which is sufficiently simple to solve. Indeed, 
the local part of the density fluctuations in a driven lat- 
tice gas is strongly influenced by an externally driven 
force [2|. However, even in this case, it has been nu- 
merically shown that the local part can be described by 
an extended, operationally constructed thermodynamic 
function [2j • This result has been proved mathematically 
for some systems under a special condition 

In this paper, in order to check the generality of the 
result in Refs. and @, we study a model of interact- 
ing Brownian particles under non-equilibrium conditions. 



Here, note that the motion of Brownian particles is be- 
lieved to be described by Langevin equations. This has 
been confirmed experimentally by measuring the trajec- 
tories of Brownian particles with a high accuracy, which 
might have become possible by the recent development 
of technology for optical instruments [f| 0, HI • From this 
fact, the system of Brownian particles is regarded as an 
ideal system for studying the fundamental problems of 
statistical mechanics, both theoretically and experimen- 
tally. 

This paper is organized as follows. In Sec. HH we intro- 
duce Langevin equations describing the motion of many 
Brownian particles under an external driving force. We 
also define the correlation coefficient C and the response 
coefficient R that we study. After reviewing briefly, the 
equality between C and R (called a fluctuation-response 
relation), we conjecture in Sec. IHII that the fluctuation 
response relation can be extended to non-equilibrium sys- 
tems on the basis of the consideration of Einstein's for- 
mula on macroscopic fluctuations. In Sec. IIV1 we report 
a result of numerical experiments. The final section is 
devoted to a few remarks. 



II. MODEL 

We study the system that consists of N Brownian par- 
ticles suspended in a two dimensional solvent of temper- 
ature T. Let Xi, i = 1,2, ...,N, be the position of the 
i-th particle, where X\ G [0, L] x [0, L] with a periodic 
boundary conditions. We express the a-th component of 
Xi as Xi a with a = 1,2. That is, X{ = {xn,Xi2). Each 
particle is driven by an external force fe\ = (/, 0) and 
is subject to a periodic potential U{x\) with period I, 
For simplicity, we assume that the periodic potential is 
independent of X2- Furthermore, we express the interac- 
tion between the i-th and j-th particles by an interaction 
potential u(\xi — Xj\). 
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The motion of the i-th Brownian particle is assumed 
to be described by a Langevin equation 



dt 



f 



dU{x z 



N 
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Sal 



u(\Xi - Xj\) + Ri a (t), (1) 



where 7 is a friction constant and Ri a {t) is zero- mean 
Gaussian white noise that satisfies 



(R ia (t)R^(t')) = 2-yT6ij6{t - t'). 



(2) 



Here, the Boltzmann constant is set to unity. The i- 
th and j-th particles interact via the soft core repulsive 
potential u(r) given by 



u{r) 



K{r - 2r c ) 2 /2 when r < 2r c , 







when r > 2r c . 



(3) 



where r c is the cutoff radius of the potential. We also 
assume a form of the periodic potential U{x\) as 



U( Xl 



Uq sin 



2nxi 

£ 



(4) 



It should be noted that without the periodic poten- 
tial U(xi), the system is equivalent to an equilibrium 
system in a moving frame with velocity f /j. Thus, the 
periodic potential is necessary for investigating the non- 
equilibrium nature. In the case that / — 0, the Langevin 
equation given in Eq. ([1]) satisfies the detailed balance 
condition with respect to the canonical distribution; con- 
sequently, the stationary distribution is canonical. On 
the contrary, in the non-equilibrium cases where / 7^ 0, 
the stationary distribution is not the canonical distribu- 
tion because of the lack of the detailed balance condition. 

In this paper, all quantities are converted to dimension- 
less forms by setting 7, £, and T to unity. The values of 
the parameters in our model are chosen as follows. First, 
r c — 2.0, that is, the interaction range of the particle is 
comparable in size to the period of the periodic poten- 
tial. Second, in the soft-core repulsive interaction u(r) in 
Eq. (HJ), K = 200,300,400, and 500. Third, U = 15.0 
and 20.0, and finally L — 30 and N = 50. Furthermore, 
Eq. (TTJ) is solved numerically by using an explicit time 
integration method with the time step At = 1.0 x 10~ 3 . 

With these parameters, wc first perform preliminary 
measurements so as to check how far the system is from 
the equilibrium. Concretely, we measure the current J 
as a function of an external force /, where we calculate 
J numerically as 



*'"+^ A dx tl (t) 



dt 



dt 



(5) 



Here, we choose the values of t ln and r c so that the right- 
hand side becomes independent of these values within the 
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FIG. 1: Current J in the xi direction as a function of /. 
The average values of J for 10 samples are plotted for several 
values of /. The statistical error bars are smaller than the 
symbols. Inset: Close-up of the linear response regime. 



numerical accuracy we imposed. The result is displayed 
in Fig. Q] by which we judge that the system with / = 20 
is in a state far from the linear-response regime. Thus, 
in the argument below, we investigate the system with 
/ = 20 as an example of non-equilibrium systems. 

Now, we study the statistical properties of the density 
field 



JV 



P(x,t) = ^2$(X - Xi(t)). 



We consider the Fourier series of the density field 



P(n 1 ,n 2 ) 



d 2 xp(x)e- i ^- {nixl+n2X2) . 



(6) 



(7) 



Using this, we characterize density fluctuations on large 
scales in terms of the correlation coefficient 



C = ((Repp (0 ,i)]) J 



(8) 



where 5pi nitn2 ) is a Fourier series of the deviation 
Sp(x) = p(x) — (p(x)) from the average density profile 
(p(x)), and Re[ ] represents the real part of a complex 
number. 

According to the fluctuation-response relation, which 
is valid in equilibrium cases, C is related to a response co- 
efficient against a potential perturbation. In the present 
problem, we consider the response of Re[<5/5( ,i)] to the 
potential 



V(x) 



e cos 



/ 27ra;2 



V L 



(9) 



We note that this potential does not depend on X\, that 
is, we apply the perturbation potential transversal to the 
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external field / in a manner similar to that in Refs. |2j,|3j. 
Quantitatively, the response coefficient is defined as 

E ^_ i , m (Re[^ (0 , 1)] ) 

e^O e 

where ( ) represents the statistical average in the steady 
state under the potential ([9|). 

Here, by the fluctuation-response relation is meant 

C = TR. (11) 

Although it should hold in equilibrium cases, there is no 
reason to expect that it also holds in non-equilibrium 
systems. Nevertheless, we can measure both C and R 
even for non-equilibrium systems. Thus, we can check 
the relation concretely by numerical experiments. 

III. CONJECTURE 

Before presenting the results of numerical experiments, 
we review the understanding of the relation given in Eq. 
(fTTj) . In equilibrium cases, the relation can be directly 
obtained from the canonical distribution; however, here 
we present an alternative, phenomcnological understand- 
ing based on a macroscopic fluctuation theory in order 
to consider an extension of the relation. 

Let p(x) be a fluctuating density field defined at a 
macroscopic scale. Then, due to the large deviation prop- 
erty, the stationary distribution is written as 

P s (p(.))^exp(-L 2 /(p(.))), (12) 

where I(p(-)) is called a large deviation functional. Ac- 
cording to the fluctuation theory for systems under equi- 
librium conditions, I is determined by thermodynamics. 
Concretely, using a free energy density f(T, p), we can 
write 

7(p(.)) = ^ / d 2 x[f(T,p(x))-f(T,p) 

+ V(x)(p(x)-p)}, (13) 

where p is the average density. This relation is called 
Einstein's formula and can be proved within a framework 
of equilibrium statistical mechanics. Noting that P(o,i) is 
an example of macroscopic density fluctuations, we find 
that Eq. (HI]) is obtained from Eqs. (Q2J) and (fT3]) . 

Equations (fT2")) and (fTB"]) are so simple that they are ex- 
pected to be extended to those valid in non-equilibrium 
systems. More concretely, we conjecture that Eq. (fT5|) 
holds if the free energy is extended in a consistent man- 
ner. 

In this conjecture, the existence of the extended free 
energy is highly nontrivial. Indeed, severe restric- 
tions are required to construct thermodynamics in non- 
equilibrium steady states, as discussed in Ref . 0] . Never- 
theless, we already know that thermodynamics can be ex- 
tended in non-equilibrium states for driven lattice gases 



when we focus on the transversal direction to the driving 
force Since the present model has common features 

with driven lattice gases, we expect that the thermody- 
namics can be extended to our model in a manner similar 
to that for driven lattice gases. 

These considerations lead us to conjecture the validity 
of Eq. (fTTj) . although the statistical distribution is not 
canonical in non-equilibrium cases. Note that Eq. (fTTj) 
can be checked independently of the question whether or 
not thermodynamic functions can be extended to non- 
equilibrium steady states. 

IV. MEASUREMENT AND RESULT 

Before measuring C and R, we estimate the relaxation 
time td for the density field. This estimation is necessary 
for the accurate determination of statistical averages by 
numerical experiments. Since perturbation effects of the 
potential V(x) can be neglected, we estimate td only for 
the systems without this potential. Initially, we prepare 
an inhomogeneous particle configuration; as an exam- 
ple, all particles are placed on the sites of a hexagonal 
lattice with the lattice interval 2r c . Then, we observe 
a diffusion phenomenon of the density field by solving 
Eq. UJ numerically. Because we found that the diffusion 
in the x\ direction is slower than that in the xi direc- 
tion, we measure only the time dependence of /5( 10 )(i). 
Repeating this procedure 30 times, we estimate the sta- 
tistical average of \f)(i,o) (t)\- We then obtain a fitting 
|p ( i )0 )(t)| = |/5 ( i i0) (0)|e-*/ rD . We found that the K de- 
pendence of td is small and that td is a decreasing func- 
tion of / when the other parameters are fixed. From 
these observations, we estimate td for the system with 
Uq, K, and / by measuring td for the system with / = 0, 
K = 200, and V(x) = 0. 

Using td evaluated above, we estimate a statistical av- 
erage {A) e using a time average 

(A) t = - f° +T dtA(t), (14) 

where we chose the values of to and r as those more than 
40 times of the relaxation time td for the density field. 

We now measure R in Eq. (fT0|) for all the values of Uq 
and K (mentioned previously) both under the equilib- 
rium condition / = and the non-equilibrium condition 
/ = 20. By measuring (Re[<5/3(Q,i)]) for several values of 
e, we estimate the regime where it linearly depends on e. 
The slope in this region yields R. As one example, in Fig. 
HI (Re[5p (0) i)]) e is displayed for U = 15 and K = 200: 
R is evaluated as 5.51 ±0.01 for the non-equilibrium case 
(/ = 20), while it is 4.64 ± 0.03 for the equilibrium case. 

Next, we measure C = ((Re[<5/S(o.i)]) 2 )o for the same 
values of the parameters. We find that for / = 20, the 
value of C shows clear deviation from the equilibrium 
value. This suggests that the stationary distribution of p 
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FIG. 3: C vs TR for several values of Uo, K, and /. Two same 
symbols correspond to the data obtained by using the same 
Uo and K for the equilibrium and non-equilibrium conditions. 



in the non-equilibrium case is not close to the equilibrium 
distribution. 

From the two independent measurements of R and C, 
we can check the fluctuation-response relation C = TR. 
The result is summarized in Fig [3J This suggests that 
the relation given in Eq. (fTTj) holds even in the non- 
equilibrium steady states to a similar extent as that in 
the equilibrium states. 

Here, we present a remark on the systematic deviation 
from the equality in Fig [3] Let us recall that the validity 
of the equality given in Eq. (|lip can be proved math- 



ematically in equilibrium cases. Therefore, each point 
in Fig. [3] should exist on the solid line within the sta- 
tistical error bar at least for the equilibrium cases. We 
conjecture that the slight but systematic deviation from 
the solid line in Fig. [3] is caused by numerical inaccu- 
racy associated with a choice of the value of At. In order 
to check it further, we measured how the deviation de- 
pends on the time step At in our numerical calculation. 
As shown in Fig [H C/(RT) approaches 1 in the limit 
At — > 0. Thus, we expect that in both equilibrium and 
non-equilibrium cases, the fluctuation-response relation 
holds with a higher accuracy if we perform numerical ex- 
periments with a smaller At. 



V. CONCLUDING REMARKS 

Before concluding this paper, we present two remarks. 
First, let us notice that the system we study is a typ- 
ical example of the so-called driven diffusive system 
[ToL [Til . It has been believed that in such a sys- 

tem the equal-time spatial correlation function generally 
exhibits a power-law decay of the type r~ d for a large dis- 
tance r, where d is the spatial dimension of the system 
[T3L [lI ] . This power-low decay is called the long-range 
correlation. In fact, it has been proved for the model 
given in Eq. ([l} that this type of long-range correlation 
appears when the interaction length between two parti- 
cles is sufficiently larger than the period £ [HI]. However, 
such a nonlocal behavior of density fluctuations was not 
observed in our numerical experiments. A plausible ex- 
planation for this apparent inconsistency might be that 
the system size we study is much smaller than that at 
which the long-range correlation can be observed, though 
the system size is so large that thermodynamic fluctua- 
tions can be argued. Further studies will be necessary in 
order to gain a clearer understanding. 

As the second remark, we address an example of lab- 
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FIG. 5: Schematic figure of a driven diffusive system under 
periodic boundary conditions. 

oratory experimental systems related to our study and 
consider the possibility of observing our simulation re- 
sults experimentally. First, since a periodic potential 
with a period of 6 /xm can be designed by using an optical 
instrument [l6[, let us assume that 1 = 1 corresponds to 
6/im. Then r c = 2.0 and L = 30 correspond to 12.0 /zm 
and 180 /im, respectively. Here, the cutoff radius r c may 
be interpreted as an interaction range such as the Debye 
screening length of the screened Coulomb potential. Ac- 
cording to Ref. [Tt} , such a long screening length can be 



realized experimentally. Note that the core radius of the 
particle is of the order of micrometers. Next, in order to 
realize the periodic boundary conditions in the direction 
of the external force, it might be a good method to place 
all the particles under a rotating optical tweezer with a 
velocity v |9j. Indeed, by observing the system in the 
moving frame v, we can confirm that the spatially homo- 
geneous force / = 7« is produced with periodic boundary 
conditions. (See Fig. 0) We expect that such an imple- 
mentation is possible due to the recent development of 
the optical technology 

In conclusion, we have demonstrated that the 
fluctuation-response relation given in Eq. (fTTj) is plausi- 
ble for many Brownian particles under an external driv- 
ing force when we focus on a direction transversal to the 
driving force. If this relation is valid for any average 
density, one can obtain a formula that relates the inten- 
sity of density fluctuations in the transversal direction 
to the derivative of a chemical potential with respect to 
the density 0. Then, by measuring the work required 
to change the system size in the transversal direction, 
one may confirm the Maxwell relation, which ensures the 
existence of a thermodynamic function extended to non- 
equilibrium steady states [1,0]. In this manner, we will 
confirm Einstein's formula for the system we study. This 
provides a realistic and nontrivial example for a frame- 
work of steady-state thermodynamics. 

We thank Masaki Sano, Yoshihiro Murayama, and Ku- 
miko Hayashi for their helpful suggestions. This work was 
supported by a grant (No. 19540394) from the Ministry 
of Education, Science, Sports and Culture of Japan. 



[*] Electronic address: 

takenobu.nakamura@aist.go.jp 

[1] R. Kubo, M. Toda, and N. Hashitsume, Statisti- 
cal Physics II: Nonequilibrium Statistical Mechanics 
(Springer, Berlin, 1991). 

[2] K. Hayashi and S.I. Sasa, Phys. Rev. E 68, 035104(R), 
(2003). 

[3] S. Sasa and H. Tasaki, J. Stat. Phys. 125, 125, (2006). 

[4] B. Derrida, J. L. Lebowitz, and E. R. Speer, Phys. Rev. 
Lett. 87, 150601, (2001); B. Derrida, J. L. Lebowitz, and 
E. R. Speer, J. Stat. Phys. 107, 599, (2002). 

[5] B. Derrida, J. L. Lebowitz, and E. R. Speer, Phys. Rev. 
Lett. 89, 030601, (2002); B. Derrida, J. L. Lebowitz, and 
E. R. Speer, J. Stat. Phys. 110, 775, (2003). 

[6] P. T. Korda, M. B. Taylor, and D. G. Grier, Phys. Rev. 
Lett. 89, 128301, (2002). 

[7] J. C. Crocker and D. G. Grier, Phys. Rev. Lett. 73, 352, 
(1994). 

[8] P. T. Korda, G. C. Spalding, and D. G. Grier, Phys. Rev. 



B 66, 024504, (2002). 
[9] V. Blickle, T. Speck, C. Lutz, U. Seifert, and C. 

Bechinger, Phys. Rev. Lett. 98, 210601, (2007). 
[10] B. Schimttmann and R. K. P. Zia, Statistical mechanics 

of driven diffusive systems (Academic Press, 1995) . 
[11] G. L. Eyink, J. L. Lebowitz, and H. Spohn, J. Stat. Phys. 

83, 385, (1996). 
[12] P. L. Garrido, J. L. Lebowitz, C. Maes, and H. Spohn, 

Phys. Rev. A 42, 1954, (1990). 
[13] G. Grinstein, D. H. Lee and S. Sachdev, Phys. Rev. Lett. 

64, 1927, (1990). 
[14] J. R. Dorfman, T. R. Kirkpatrick, and J. V. Sengers, 

Ann. Rev. Phys. Chem. 45, 213, (1994). 
[15] T. Nakamura and S.I. Sasa, Phys. Rev. E 74, 031105, 

(2006). 

[16] C. Bechinger, M. Brunner, and P. Leiderer, Phys. Rev. 

Lett. 86, 930, (2001). 
[17] A. Yethiraj and A. v. Blaaderen, Nature 421, 513, 

(2003). 



